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Abstract 

The full 0{as) transverse momentum distribution for the b sj decay is computed. 
Results are presented in analytic form. An improved expression for the coefficient function 
taking into account subleading operators is given and an exact expression for the remainder 
function associated with the leading operator Oy is also derived. 
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1 Introduction 



Perturbation theory, i.e. the expansion in powers of as, has been appHed to describe decays of the 
beauty quark since its discovery. While the expansion parameter as{rnB) ^ 0.21, being reasonably 
small, allows one to have confidence in the computations, it is difficult to directly compare the per- 
turbative approach with the experimental data. As is well known, decay rates do not make good 
quantities to be compared with the data, because they are proportional to the fifth power of the 
beauty quark mass, a poorly known parameter, 

r oc (1) 

and because they involve in principle unknown CKM matrix elements such as Vet, Vub, Vts, etc.. By 
taking ratios of different widths, one can cancel the dependence in the observables, and, eventually, 
also the dependence on the CKM matrix elements. A rather good theoretical quantity is represented, 
for instance, by the semileptonic branching ratio: 

BsL^^, (2) 

i TOT 

which turns out to be marginally in agreement with present data Inclusive quantities, Bindusive, 
such as have a perturbative series that involves numerical coefficients c„ of the form: 



E 



c„ asims)- (3) 



In less inclusive quantities, additional dynamical effects appear, due to the kinematical restrictions on 
the final particles, and the use of perturbation theory is, in general, less justified. In semi-inclusive 
quantities, B semi-inclusive , such as threshold and transverse momentum pt distributions, the pertur- 
bative series contains large infrared logarithms in addition to the coefficients c„; they may be expanded 
as a perturbative series of the form: 

oo 2n 

B semi— inclusive ^ ^ ^ ^ Cn,k (4) 

n=0 k=Q 

where x represents the characteristic scale of the process as the energy or the transverse momentum. 
Resummation of such enhanced terms to any order in can be performed in various approximations. 
The simplest one, the leading logarithmic approximation, involves picking up only the terms having 
two powers of the logarithm for each power of the coupling, i.e. k = 2n. In the double-logarithmic 
approximation each parton is dressed with a cloud of soft and collinear gluons. Further, more refined 
approximations involve smaller numbers of logarithms for each power of as, i.e. /c = 2n— 1, 2n — 2, . . .. 
In the last years, considerable effort has been devoted to the study of various spectra in B decays in 
the endpoint region, in the framework of resummed perturbation theory. 

In order to verify the ability of the resummed perturbation theory to describe B decays in a different 
dynamical situation, we considered, in a previous note 2j, pt-distributions describing that of the s 
quark with respect to the photon direction, in the b rest frame. 

In this work, 2 , the following issues have been considered: the resummed pt-distribution in the 6 — s- 57 
decay is evaluated and both perturbative and non-perturbative sources of transverse momentum con- 
tributions discussed. The general theoretical framework for the evaluation of the corresponding matrix 
element defined and the strategy to evaluate leading and next-to-leading perturbative contributions 
is outlined, by introducing a method to treat the radiative corrections and their summation in a im- 
proved perturbative formula. The comparison of the transverse momentum distribution singularity 
structure with the more widely-known threshold case is also presented. 

The chosen quantity manifests a clear advantage from a phenomenological point of view since, as dis- 
cussed in PI, it depends only on the photon momentum in the process. Thanks to the straightforward 
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and direct kinematics, the transverse momentum turns out to be a particularly simple variable to use 
to discuss the singularity structure of the perturbative expansion. The case of a possible effective 
theory within which to factorize these singularities can, for the transverse momentum, be considered 
as well. 

The general formula representing the complete perturbative expression for a the resummed distribu- 
tion is given by the formula 

D{x)^K{as)nx\as)+R{x\as)- (5) 

The results, already presented in [2], did concern the universal process-independent function S(a;; as), 
resumming the infrared logarithms in exponentiated form. 

Here the general perturbative expression for the whole distribution will be concisely recalled and the 
new entries represented by the coefficient function K(as) and by the remainder function R(x) will be 
evaluated. Both K(as) and R{x) are process-dependent and require an explicit evaluation of Feynman 
diagrams. 

Resummation of large infrared logarithms in b decays have been studied in great detail in recent years. 
This scheme is justified by the fact that the double logarithm appearing to order as can become rather 
large (with respect to 1 coming from the tree level): 

-^10^4 -'0.7 (6) 

if we push the transverse momentum to such small values as pt ~ ^qcd — 300 MeV. The single 
logarithm can also become rather large, having a large numerical coefficient: 

log 4 -0.6. (7) 

The purpose of resumming classes of such terms therefore seems quite justified. If we consider running 
coupling effects, i.e. if the (frozen) coupling evaluated at the hard scale Q ~ ms — 5.2 GeV is 
replaced by the coupling evaluated at the gluon transverse momentum, 

Qs(mb) asipt) = 0.45 for pt = 1 GeV, (8) 

the logarithmic terms have sizes of order: 

-^^(^log^4--0.5 (9) 
47r 

and 

5asCF . pi ... 
log— -0.8. 10 

An 

The main difference with respect to resummation in decays is a hard scale smaller by over an order 
of magnitude, i.e. a coupling larger by a factor 2 and infrared logarithms smaller by a factor 3. 



2 The effective hamiltonian for the decay 6 ^ 57 

The decay 6 ^ 57 is loop-mediated in the Standard Model and offers stringent tests of the latter as 
well as a way to extract CKM matrix elements. The relevant diagrams involve a loop with a virtual 
W and an up- type quark (u, c or i); the external photon can be emitted from the internal lines and 
from the external lines of the 5 or s quark (see fig. QCD radiative corrections are affected by large 
logarithms of the form 

log'' ^ with < fc < n (11) 
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Figure 1 : Vertex for & ^ 57 in the Standard Model 



as well as logarithms of mt/mw. Since the energies involved in the process are much smaller than the 
W OT t mass, it is possible to integrate out these fields by means of an operator product expansion 
and write an effective low-energy hamiltonian of the form: 

-Heffix) = % VZVtt V C.ifit) (12) 



With a factorization scale fib = 0{mb), the long-distance effects — both perturbative and non- 
perturbative — are factorized in the matrix elements of the operators Oj, while the short-distance 
effects are contained in the coefficient functions Cj{fib), calculable in perturbation theory. In particu- 
lar, the large logarithms in 1)11(1 are included into the coefficient functions and can be resummcd with 
standard renormalization group techniques. 

A suitable basis for the operators Oj is given by six four-quark operators, Oi-Og and by the penguin 
operators 67, Os HE]: 

Cl = ^L,plpibL,a){sL,alpiCL,l3) 
O2 = {cL,a"ftibL,a)isL,l3Jtj.CL,l3) 

9 

q 

O5 = {sL,al,,bL,a)C^qR,f3l^LqR,l3) 
9 

Oe = {sL.alf.bL.^pjC^qR^fjlt.qR.a) 
Q 

6 

O7 = Y^mi^j^{fib)sL.a(7^''bR^aFfj,i, 

^8 = Y^'^b.'MS^^^I')sL.aCr^''T^i3bR,aG1;^^, 



(13) 



where fnj^jjg{fib) is the b mass in the MS scheme, evaluated at fib and q = u,d,s,c or b. 

The dimension of these operators is six: higher-dimension operators have coefficients suppressed by 

inverse powers of the masses of the integrated particles {t and W) and do not contribute in first 
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approximation. 

The calculation of the QCD corrections to the coefficients functions has been carried out in ^ with 
leading logarithmic accuracy and in |7] at next-to-leading level in the MS scheme. 
Let us now consider the evaluation of the matrix elements of the effective hamiltonian between quark 
states. Only the magnetic penguin operator contributes in lowest order with a rate:^ 

To ^ 32^-^ C, {^^,) , (14) 

where m;, is the pole mass of the b quark. 

Radiative QCD corrections involve gluon brehmsstrahlung. The operator is affected by infrared 
singularities for the emission of a soft or a coUinear gluon; the remaining operators Oi-Oq have infrared- 
finite matrix elements. This implies that QCD corrections to the operator O7 only are logarithmically 
enhanced for pt ^ ^. We will then consider at first only the operator O7. 



3 Transverse momentum distribution in 6 ^ S7 

The process we are dealing with has a very simple kinematics: in lowest order it is the two-body decay 
b ^ sj. Let us define 

x=^, (15) 

where pt is the transverse momentum of the strange quark with respect to the photon direction, fixed 
as z-axis, and mi, is the mass of the heavy quark, to be identified with the hard scale of the process'^. 
In lowest order the transverse momentum distribution then is 

— ^ToS{x), (16) 

that is the strange quark and the photon are emitted in opposite directions, because of momentum 
conservation. AcoUinearity is generated by gluon emission; in 6 — s- sjg, i.e. at 0{as), Pt = —kt while 
in & s^gi ■ ■ ■ gn, i-e. in higher orders, pt = —kti ... — ktn- 

Beside the differential distribution the partially integrated distribution^ is also of interest 

Even though asims) is small enough to justify a perturbative approach, the combination ag{mB) log*^ x 
with < k < 2n can be large. A resummation, to any order in as, of logarithms of the same magni- 
tude is required to obtain sensible physical results. 

A partial resummation of large logarithms with next-to- leading accuracy has been performed in [5]: 
here we complete the calculation. 

It is well known, |12 | , that the resummation of large logarithms is accomplished by an expression of 
the form: 

D{x) ^ K{as)^ix;as) + R{x;as), (18) 

where 



^Fq contains in principle m'^——{iii,) since ni, is the renormalization point of O7. As is well known, the renormal- 
ization point is arbitrary: we decided to fix it to mj, in the running mass, as usually done in the literature. 
■^The operator Og is affected by QED infrared divergences which are not relevant to our problem. 
^Let us note that < a; < 1/4. 

^Since we divide the spectrum by the lowest-order rate Fq, we have that D(x = 1/4) = ^TpT = 1 ^ 0{ag). 
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• 'E{x;as) is a universal, process-independent, function resuniming the infrared logarithms in 
exponentiated form. It can be expanded in a series of functions as: 

logT,{x;as) ^ Lgi{asL) + g2{asL) + asgsiasL) + . . . , (19) 

where L = logx (in general L is a large infrared logarithm). The functions gi have a power 
expansion of the form 

oo 

g,(z) = ^g,,fc z'^ (20) 

fc=0 

and resum logarithms of the same size: in particular gi resums leading logarithms of the form 
Og L""*"^ and g2 the next-to-leading ones Og L". The explicit form of S(a;; as) can be found in 
Ref.0; 

• K{as) is a short-distance coefficient function, a process-dependent function, which can be cal- 
culated in perturbation theory: 

if (as) = 1 + ^^fci -I- O (4) . (21) 

• R{x]as) is the remainder function and satisfies the condition 

R{x] as) for x 0. (22) 

It is process dependent, takes into account hard contributions and is calculable as an ordinary 
as expansion: 

R{x-as)^^^r^{x)+0{al). (23) 

The result is an improved perturbative distribution, reliable in the semi-inclusive region T?, that is 
for small values of x, which can be matched with a fixed-order spectrum, describing the distribution 
for large values of a; 13 . The description of the tools used to perform the rcsummation of infrared 
logarithms is far from the purpose of this note, and we refer the reader to the references |S] - |12| . 

In the next sections the full order as corrections for the coefficient function K{as) and for the 
remainder function R(x;as) will be explicitely calculated. 



4 0{as) corrections to pt distribution 

In this section radiative corrections to the transverse momentum distribution will be calculated eva- 
luating the Feynman diagrams depicted in fig. |21 for real gluon emissions and in fig. |31 for virtual 
emissions. We use the Feynman gauge where the gluon propagator is 

D,Ak) = (24) 
The calculation is performed in dimensional regularization (DR) with the dimension of the space-time 

n = 4 + e. 

The operator O7 from the basis (|13|l is inserted in the hard vertex, as discussed in section|21 Let us now 
define the kinematical variables: is the heavy quark momentum, p'' the light quark momentum, fc^ 
the gluon momentum and the photon momentum: for real diagrams it holds that = = = Q 
and — m^, while, for virtual diagrams, fc^ ^ 0. The calculation is performed in the b rest frame, 
where 

= (i?^,0,0,£;^). (25) 
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Figure 2: Real diagrams 



4.1 Real diagrams 

A straightforward evaluation of the diagrams in fig. |3 gives a contribution to the rate: 



dT _ M{uj,t;e) 



To oj^-H^-'/^ 



dt dco — 



OJ 



l-£^l-c/2 



dt duj 



(26) 



where 



Si it) 



and 



= M(0,0;e) 

M(0,i;e) - M(0,0;e) 

= 

A/(w,0;e) -M(0,0;e) 

Mjuj, t- e) - A/(0, t; e) - M{uj, 0; e) + M(0, 0; e) 
^ 1 — cos 9 



(27) 



2F • fc 2E„ 



rrib 



with the angle between the gluon and the direction —z^. 

It follows from their definition that the functions Ai(w, t\ e), S'i(w; e), Ci(i; e) and -Fi(w, e) are finite 
in the soft and the coUinear limit, defined respectively as 

cj and 0. (28) 

The rate in eq. H27|l has to be integrated over the whole phase space with the kinematical constraint 



5[x ~ Loh{l ~ t)], 



(29) 



^Let us remember that the direction +z is fixed by the photon space momentum. 
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which selects gluons with transverse momentum = xml ^. The cumulative distribution takes a 
contribution of the form: 



/■! /■! I dY 

Dr{x)= i dx' I dhj I dt — —{Lo,t-e) 5[x' -Loh{l-t)]. (30) 
Jo J a Jo ^0 dx 

After the integration we expect four kinds of terms: 

• Poles in the regulator e: they parametrize the infrared singularities and cancel in the sum with 
virtual diagrams because the distribution we are dealing with is infrared-safe''; 

• Logarithmic terms diverging for a; — > 0; 

• Constant terms: they enter the coefficient function K{as)] 

• Remainder functions: terms that vanish in the limit a: 0. 
Integrating over x' we have: 

1 '•1 1 dV 



Dr{x)= du dt — —{uj,t-t) e[x-ujH{l-t)]. (31 



Jo To dx 

The remaining integrations are non-trivial because of the simultaneous presence of the kinematical 
constraint and by the dimensional regularization parameter e. By using the identity 

e[x - ujh{i - 1)] = i-e[uj^t{i-t)-x\, (32) 

we separate these two effects and rewrite the distribution Dr(x) as a difference between an integral 
over the whole phase space and a integral over the complementary region: 



/•I 1 dV I dV 

Dr{x)^ duj dt— —{uj,t-e) - / du dt— —{u;,t;0) 0[Loh{l 
Jo Jo To dx Jo Jo To dx 



t)~x] + 0{e). (33) 



The first integral must be evaluated for e because it contains poles in e, but is done over a very 
simple domain, independent of x. The second integral does not contain any pole in e and therefore 
one can take the limit e ^ in the integrand. It depends on the kinematics of the process and can be 
integrated by introducing a suitable basis of harmonic polylogarithms as in |14|. The most convenient 
basis we found consists of the basic functions 



<?[-!; 2/] 
<?[-2;2/] 



y + l 

1 



The harmonic polylogarithms (HPL) of weight 1 are defined as: 



J[a;y] = / dy' g{a;y') for a 7^ 
Jo 

J[0;y] = logy. (35) 



®Let us recall that for the single gluon emission x = p^/m^ = k^/m^. 

distribution is infrared-safe if it is insensitive to the emission of a soft and a coUinear gluon 1131 . 
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In terms of usual functions, they read: 



J[-l;y] = log(l + y) 
J[-~2-,y] = 2arctany^ 
•^[-3;?;] = log(l-V^7^) 

The HPL's of weight 2 are defined for (u, v) ^ (0, 0) as 



(36) 



J[u,v;y]= dy' g[u;y'] dy"g[v;y" 
Jo Jo 



(37) 



and J[0, 0; y] = 1/2 log^ y. HPL s of higher weight may be defined in an analogous way. They will not 
be used here. 

The final result for real diagrams turns out to be 



as 



e/2 



1 



r(I + e/2) 



5_ 
Ye 



1 5 1 

- log^x- - logxH- - +d{x) 



(38) 



where d{x) is a function vanishing for a; — > 0. 

The matrix elements of the remaining operators Oi (i ^ 7) do not contain infrared divergences. 
Therefore their contributions to Dji do not involve (infrared) poles in e, logarithms of x and constants, 
but only new functions, which vanish in the limit x ~* 0. 

4.2 Virtual diagrams 

Virtual corrections to b —^ sj have been calculated in '15', '17' for a massive strange quark; we present 
here the computation in the massless case. The diagrams consist of self-energy corrections to the 
heavy and light lines and of vertex corrections to the operator Or (see fig. OJ; we compute them in 
the MS scheme so as to be consistent with the (known) coefficient functions Ci. The computation 
can be done with standard Feynman parameter technique or by a reduction using the integration by 
part identities ^H]- Let us briefly describe the evaluation of the vertex correction within the latter 




Figure 3: Virtual diagrams 
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method. One has to compute the scalar integrah 

f _d-k_ N{k\P-k,p-k-t) 

J (27r)" k-^[{k~ Py -m^]{k~pY' ^ ' 

where 

N{k^,P ■k,p-k;e) = 32P ■ k - 32p ■ k - 16mg + 0{e^)k^ + 0{e^)P ■ kp ■ k + 0{e^){p ■ kf. (40) 

V has at most a double pole in e coming from the product of the soft and the coUinear singularities. 
The terms in the numerator TV, which vanish in the soft limit /c^ — > 0, do not give rise to soft 
singularities and therefore produce at most a simple pole coming from the collinear or the ultraviolet 
region. Therefore the O(e^) terms in do not contribute in the limit e — > 0. 

By expressing the scalar products in the numerator as linear combinations of the denominators as 

k-p ^ ^{^-{k-pf), 

k-P - ^{k^-{k-Py + ml), (41) 

we can reduce V to a superposition of scalar integrals of the form: 

r dJ'k 1 
Tia. ^' - y [fc2]a [(^_p)2_^2p [{k-pYY 

with a,b,c < 1. The above amplitudes can be related to each other by identities of the form |18|: 

d v^^ 

f]ni. _r _ n (A'i] 

dkf^ [k^Y[{k - Py - m^]b[{k ^ p^Y 

with — k'^jp^jP'^. By explicitly evaluating the derivatives and re-expressing the scalar products 
using eqs. (|41l) . one obtains relations among amplitudes with shifted indices. 

By solving the above identities, one can reduce all the amplitudes to the tadpole, and one obtains for 
our integral*: 

V= f-— + 8-86) ^ T[0, 1,0], (44) 

where 

Summing self-energies and vertex corrections, and subtracting the 1/e poles according to the A/5 
scheme, one obtains for their contribution to the rate Dy ^■ 



(^S ml \ / e 



2 5 mb 
- + - + 41og — 
2e fib 



(46) 



We have kept the factor in front of the square bracket unexpanded to simplify the computation of the 
total rate. 

The virtual corrections to the remaining operators O^^y contain only (simple) ultraviolet poles in e, 
which are removed by renormalization; their contributions to Dy amount only to finite constants and 
logmf,/^f,. 

*Such a strong reduction of 3-point function to a vacuum amplitude is possible because the only scale in the process 
is the heavy quark mass rrii,. Virtual corrections have indeed the lowest-order kinematics = m^, P ■p = = 
(j2 = 0. 

factorize Fq one has to replace m^-^j-^{^i,) by m^j^{mi,) using the formula m^Yf§{fii,) = m^j^{rrn,){l + 
I). 



3 Cpas ' 
2 
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4.3 Final result 



Summing real and virtual contributions, the transverse momentum distribution for the decay 6 — > 37 
reads, to 0{as). 

"1,2 5 



D{x) = 1 + Cf — 

TT 



log X - - \ogx + f+ d{x) 



(47) 



As expected, the result contains a double logarithm and a single logarithm of a;, a finite term / and 

a function d{x) vanishing in the limit x — s- 0. 

By expanding the resummed formula to order as one obtains: 

Cpas , \ f, Ai 2 _ , n _ 1 \ , Cpois 



D{x) = \l + —^kij\l-—as\og''x + Bias\ogxj+—^r{x) 

= l-4ia5log'a; + Siasloga; + ^^fci + ^^r(x)+0(a|). (48) 

4 TT TT 

By identifying the resummed result expanded to 0{as) with the fixed-order one — matching procedure 
— we check the values for Ai and Bi evaluated in our previous paper using general properties of QCD 
radiation |2] and we extract the value of the coefhcicnt function: 

as well as the remainder function ri{x) — d{x). As explained in previous sections, the remaining 
operators O^^y contribute to D(x) only by finite terms fi and remainder functions. Since the constants 
fi come from virtual diagrams alone, we can quote their result from jl5| and present an improved 
formula for the coefficient function, in analogy with |16| : 

where 

fi ^ ri i y^7 

rr = ^f/-41og^U-^-^. (51) 

Let us remark that only the coefficients related to the operators with i = 1, 2, 7, 8 are relevant, because 
the others are multiplied by very small coefficient functions and can be neglected: 

1 

n = --r2 
6 

sRra = -4.092- 12.78(0.29 -TOc/mf,) 

rs = ^{33-2tt'). (52) 

The analytic expressions for the coefficient functions as well as a standard numerical evalutation are 
given in 7 . The anomalous dimension "f!ff is derived from the coefficient of the logarithmic term in 
ki . The values of 7^7 ^ are [3 : 



(0) _ / 208 416 176 152 6272 4624 32 32', 



Eauation()50|l is the main result of our paper and allows a complete resumniation to NLO of transverse 
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momentum logarithms. 

The expHcit calculation of the remainder function in 1)48(1 reads 



log 



_ (t - l)(49r* + 468t'^ + 1797r6 + 3642t^ + UbOr^ + 3642r3 + ITOTr^ + 468r + 49) 
" 12(r + l)5(T2+3r + l)2 

-5 - 61t - 317r2 - 912t^ - IQ22t^ - 1934r5 - IQ22t^ - 912r^ - 317t« - 61r^ - Sr^o 

4(t + 1)6(t2+3t + 1)2 

J[0, -3, t] + J[0, -3, l/r] - 2 J[0, -1, t] + J[-l, 0, t] + J[-l, -3, t] - J[-l, -3, 1/r] 

^ _ , „ (t + 1)(2t2 + 7t + 2) tt _(r + l)(2T2 + 7r + 2) 49 

2V7arctan(V7)^-pL__^ + S^^^TsTTT)^ + 

^logr- ^log(T + l)+log2(r + l), (54) 



where 



r=l^^^. (55) 

1 + Vi^i 

Let us notice that r behaves as x for small values of the transverse momentum 

t{x)=x + 0{x^) (56) 



and it is a unitary variable 



T ^ for a: ^ 
T ^ 1 for a: ^ 1/4. 



The relation 1(55(1 may be inverted as 



(57) 



(r + l)2- 

One can easily check that r(T) vanishes for t x by using the properties 

J[0, -1, 0] = J[0, -3, 0] = J[-l, 0, 0] = J[-l, -3, 0] = (58) 



7r2 



lim J[0,-3,l/-r] = lim J[-l,-3, 1/t] = - — . (59) 

T— »0 r^O 3 



5 Conclusions 



To sum up: eq. l|T7)l . which contains the final result, represents the full evaluation to 0{as) of the 
transverse momentum distribution. It is explicitly given in terms of an analytic expression. 
Contrary to what happens in hard processes at much larger energies, at the energy scales involved 
here for the h decay the remainder function contribution does play a more important role. 
A straightforward numerical evaluation of the remainder function r{x) of eq. I(48() allows us to conclude 
that its contribution can be safely neglected for small values of x, up to x ~ 0.1, where it approaches 
the zero limit of a; ^ 0. For larger values of a:, however, the size of its contribution increases to reach 
values of the order of the 10-15% of the combined leading and next-to-leading logarithmic terms. 
A detailed report describing the calculation giving rise the 0{as) evaluation presented here, together 
with an analysis of the related phcnomenological impact, will be presented in a future article [3]. 
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